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Abstract. A new Hamiltonian structure of the Maxwell-Bloch equations is described. In 
this setting the Maxwell-Bloch equations appear as a member of a family of generalized 
Maxwell-Bloch systems. The family is parameterized by compact semi-simple Lie groups, 
the original Maxwell-Bloch system being the member corresponding to SU(2). The Hamil- 
tonian structure is then used in the construction of a new family of symmetries and the 
associated conserved quantities of the Maxwell Bloch equations. 
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1 Introduction 

Maxwell-Bloch equations are a system of partial differential equations which plays a prominent 
role in the field of non-linear optics. This system models the resonant interaction between 
light and an optically active medium consisting of two-level atoms. The quantities figuring 
in this system are the complex valued functions E(t,x), P(t,x): R 2 — > C and a real valued 
function D(t,x): R 2 — > R of two independent variables, time t and one spatial variable x. The 
function E(t, x) is the slowly varying envelope of the electric field, P(t, x) is the polarization 
of the medium and D(t,x) is the population inversion. Here we will consider the Maxwell- 
Bloch equations without pumping and in the sharp-line limit, that is, without inhomogeneous 
broadening: 

E t + cE x = P, P t = ED-(3P, D t = -^(EP + EP). (1) 

The constant c above is the speed of light in the medium and (5 represents the longitudinal 
relaxation time of the medium. The transverse relaxation time and the losses of the electric 
field are assumed to be equal to zero. 

We will show that the Maxwell-Bloch equations describe a continuous chain of interacting 
C. Neumann oscillators on the three-sphere. The interactions between the neighbouring oscilla- 
tors are of magnetic type, which means that the acceleration of any given oscillator depends on 
the velocity of its neighbours and not on their position. More explicitly, we will show that the 
system (^Q) is equivalent (modulo a certain constraint) to the second-order partial differential 
equation 

{gtg- 1 ) t + c{g t g- 1 ) x = [a,M g (r(x))], (2) 

where g(t,x): Kxl^ SU(2) is a Lie group valued function of two variables, a £ su(2) is 
a constant and t{x) : R — > su(2) an arbitrary path in the Lie algebra. We shall see that this 
equation is the equation of motion of the above-mentioned chain of C. Neumann oscillators. 
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If we replace the group SU(2) by an arbitrary Lie group G, the equation @ still makes sense. 
We shall call the equation (J2J), with g(t, x): Ixl^G and a, t(x) G = Lie(G), the generalized 
Maxwell-Bloch equation. We will show that in the case when G is a compact semi-simple Lie 
group, the equation (J2J) is again the equation of motion for a continuous chain of oscillators. 
In this case, the configuration space of the relevant model oscillator will be the Lie group G. 
These oscillators belong to a class of well-known integrable systems described by Reyman and 
Semenov-Tian-Shansky in and and later in different contexts by other authors, see e.g. [3]. 

The rewriting (J2J of the Maxwell-Bloch equations proves to be useful in several ways. In 
this paper, we shall concentrate on two features. First we shall construct and describe a new 
Hamiltonian structure for the Maxwell-Bloch equations and their generalizations. This will 
then enable us to find an infinite-dimensional group of symmetries and corresponding infinite- 
dimensional class of conserved quantities of our equation. 

We shall limit ourselves to the spatially periodic case 

g(t, x + 2ir) = g(t, x), t(x + 2ir) = t(x), 

and we shall assume that the group G is compact and semi-simple. We will show that (J2J is the 
equation of motion for the Hamiltonian system (T*LG, oj c + cio m , H^). Here the configuration 
space LG = {g(x): S 1 — > G} is the loop group over G. The Hamiltonian H m b'- T*LG — » R is 
given by the formula 



where K(—,—) is the Killing form on g. The symplectic structure oj c + cui m is the canonical 
structure oj c on the cotangent bundle T*G, perturbed by the so called magnetic term cuj m . The 
term u> m is the pull-back 7r*(c5 m ), via the natural projection ir : T*LG — > LG, of a right-invariant 
differential 2-form ui m on LG. The value of ui m at the identity is given by 



We stress that in the case when G = SU(2), the above system provides a new Hamiltonian 
structure for the Maxwell-Bloch equations (0). 

We see that the symplectic structure of our Hamiltonian system is not canonical. It is well- 
known that the perturbations of the canonical symplectic forms are responsible for the forces 
of the magnetic type, see e.g. [UOHIE]. The model example of such a situation is the motion 
of an electrically charged particle in a magnetic field. The Lorentz force can be encoded as 
a perturbation of the canonical symplectic form on T*R 3 . One can geometrize the Lorentz-type 
force by adding another (circular) degree of freedom to the configuration space. On the suitably 
extended phase space T*M the symplectic form will be canonical. One can then easily find 
such a metric on M that the geodesies on M will project down to the trajectories of a particle 
in N under the influence of our magnetic-type force. This procedure is known by the name 
of the Kaluza-Klein theory. Whenever the magnetic perturbation uj m is exact, the Kaluza- 
Klein extended space is simply M = N x U(l). If uj m is not exact, then the extended space 
exists under the condition that u> m has a certain integrality property. In this case, the extended 
space is a non-trivial [/(l)-bundle whose Chern class is equal to the de Rham class of u) m . In 
the case of the generalized Maxwell-Bloch equations the form u m G Q 2 LG is not exact. The 
extended configuration space is therefore a non-trivial L r (l)-principal bundle over LG. It is 
actually precisely the central extension LG of the loop group LG. A detailed exposition of the 
Kaluza-Klein description of the Maxwell-Bloch equations can be found in [Jj. The geometric 
prerequisites, needed for the construction of non-trivial extended spaces, are given in [Hj. Here, 
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we would only like to mention the following interesting fact: In the Kaluza-Klein description 
of the electron moving in a magnetic field, the moment, conjugate to the additional circular 
degree of freedom, is the charge of the electron. Therefore, also in other situations, this moment 
is called charge. It is interesting to note that in the case of the Maxwell-Bloch equations the 
charge is precisely the speed of light in the medium. 

As we already mentioned, the important merit of the equation (J2J) and the corresponding 
Hamiltonian system (T*LG,u c + cw m , H m b) is the fact that they enable us to find an infinite- 
dimensional family of symmetries and the corresponding conserved quantities of the (generalized) 
Maxwell-Bloch equation. We will show in the Section 0] that an Abelian loop group LT, where 
T is a certain maximal torus in G, acts in a Hamiltonian way on our system. The corresponding 
conserved quantities will be constructed as the components of the moment map associated to 
the action of LT. Thus, we shall obtain a new family of conserved quantities. 

The rewriting © of the Maxwell-Bloch equations appeared for the first time in the papers [0] 
and jlUj by Q-Han Park and H.J. Shin. Without previous knowledge of these papers, the author 
rediscovered this rewriting and used it in the paper 0. In [§] and J0J the equation (J2J) is viewed 
as an equation of field theory, while in [7] it is treated as a continuous chain of C. Neumann 
oscillators. With respect to the generalized Maxwell-Bloch equation it is interesting to note 
the following. The authors of jlUj show that certain cases of the generalized Maxwell-Bloch 
equation (j2j) with SU(2) replaced, say, by SU(3) or 5C/(4) are actually physically meaningful. 
They describe various instances of degenerate and non-degenerate interaction of light with two 
and three level optical media. The correct choice of r(x) is crucial here. 



2 Generalized Maxwell-Bloch equation 

First we will show that the equations and © are indeed equivalent. Let us introduce the 
matrix valued functions p(t, i):Rx!-> su(2) and F(t, x) : R x R — > su(2) by the formulae 

In terms of F and p the Maxwell-Bloch equations become 

fH = [ P ,F], F t + cF x = [p,a], (4) 

where 

1 (i 
o = - 



2 \0 -i 

is the first Pauli matrix multiplied by i. The first of the above equations is a Lax equation. Its 
general solution is of the form 

p(t,x) = Ad 9 ( ti!C )(r(a;)), F(t, x) = -g t (t, x) ■ g' 1 ^, x). (5) 

Here g(t, x): MxM-> SU{2) is a Lie group valued function and t(x) : R — > su(2) takes values in 
the Lie algebra su(2). If we insert (JSJ) into the second equation of @, we indeed get the second 
order equation 

{9t9- 1 ) t + c(g t g- 1 ) x = [a,Ad g (r(x))]. (6) 

We note that the diagonal terms of the matrix F(t, x) are constant. Thus we see that Maxwell- 
Bloch equations are equivalent to the equation © together with the constraint 

(g t g~ l ,(r) = const = -/?. 
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From the point of view of the original physical interpretation, it would be better to introduce 
the variables p and F in the formula © as Hermitian rather than skew-Hermitian matrices. 
This would be achieved by multiplying the matrices by —i. Then p would indeed have the 
proper form of the projection on the quantum mechanical wave-function. We have adopted the 
skew-Hermitian form because it is more convenient from the mathematical point of view on the 
one hand, and because it is better suited to our mechanical interpretation of the Maxwell-Bloch 
equations on the other. Namely, being skew-Hermitian, the matrices p and F are the proper 
elements of the Lie algebra su(2) and this will facilitate our reasoning and calculations below. 

From now on we shall denote the Killing form by (—,—). In the equation @ we can replace 
the unknown function g(t, i):KxK-> SU(2) by the function g(t, i):lxl-»G which takes 
values in an arbitrary Lie group G, if only we replace a £ su(2) and t(x) £ su(2) by an element 
a € g and a curve t(x) : R — ► g, where g is the Lie algebra of G. 

We can simplify the equation to some extent. First we note the following. Let g(t, x) : R x 
R — > G be a solution of ©. Then it is easy to check that for every h(x) : R —> G the function 
f(t, x) = g(t, x)h(x) : R x R — > G is a solution of the equation 

{ftr l )t + c{f t r 1 ) x = [a, Ad f {h(x) ■ r(x) ■ h'^x))] . 

Therefore, we can assume that t(x) : R — > t C g is a map which takes values in a chosen maximal 
toroidal subalgebra t C g. 

We obtain an important equation, if we perform the following reduction. Consider the equa- 
tion @ for the case when G = SU(2), but let the unknown function g(t,x) be constrained to 
take values only in a subgroup U(l) C SU(2): 

- (-ltd) ::: ( ( r4) um c su ^- 

Let a be the Pauli matrix a = i diag (i, —i) and let t(x) = a. Then a calculation shows that 
the equation © gives 

<Ptt + C(f) tx = sin cj). 
This equation is essentially the sine-Gordon equation. 



3 Hamiltonian structure 

Let us consider the spatially constant solutions g(t) : R — > G of the equation @. These solutions 
actually solve the equation 

( 9t g- 1 ) t = [vM a {T)]. (7) 

The following proposition describes the Hamiltonian nature of the above ordinary differential 
equation. 

Proposition 1. The equation (J7J) is the equation of motion of the system (T*G,lo c , H), where 
u> c is the canonical symplectic form on the cotangent bundle T*G and the Hamiltonian function 
H : T*G — > R is given by 



H(g,p g ) = ^\\p g \\ 2 + (a,Ad g (r)). 



(8) 
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Proof. Let us trivialize the cotangent bundle T*G over the Lie group G by means of the right 
translations, T*G = G x g. With this trivialization in mind we shall denote the elements in the 
tangent spaces T^ gpg ^(T*G) = g x g* by (X b ,X ct ). The canonical symplectic form uj c on the 
cotangent bundle T*G over the Lie group G can be given by the formula 

(u c ) {g>Pg) ((X b ,X ct ), (Y b ,Y ct )) = -(X ct ,Y b ) + (Y ct ,X b ) + ( Pg , [X b ,Y b ]). (9) 

Above (— , — ) denotes the pairing between the elements of g and those of g*. For the proof of 
this formula see jllj . 

A solution of a Hamiltonian system is an integral curve of the Hamiltonian vector field Xh 
and thus is given by the relation dH = ui(Xh,—), where u is the symplectic form. For the 
Hamiltonian given by (jSJ we have 

{dH,(5 b ,5 ct )) = -([a,Ad g (T)]\5 b ) + (5 ct ,pl), (10) 

where b: g — > g* and (J: g* — > g are defined by a 9 = (a,—) and (3 = ((3$,—). Let us denote 
Xh = (X b , X ct ). Then we have 

{uc){g,p g ){.{X b ,X ct ),(5 b) 5 ct )) = -(X ct ,5 b ) + (5 ct ,X b ) + (p g , [X b , S b ]) 

= (-X4 ~ {X b , Pg },5 b ) + (8ct,X b ), (11) 

where {a, a} denotes the ad*-action of a S g on a S g*. The above equations now yield 

p\ = X b , [a, Ad g (r)f = X ct + {X b ,p g }, 
and therefore 

X b =pl X ct = [a,Ad g (T)] b . 

Let now 7(t) = (g(t),p g (t)) : M. — > T*G be a curve expressed in the right trivialization and let 
7 = {9t9~ l , {Pg)t) be its tangent at the point (g(t),p g (t)). The above equations tell us that ^(t) is 
an integral curve of the Hamiltonian vector field Xh, if and only if it satisfies the equation 

(gtg-^t = k,Ad fl (r)]. ■ 

In the case when G = SU(2), the system (T* SU(2),uj c , H) is the C. Neumann system which 
describes the motion of a particle on S 3 = SU (2) under the influence of a quadratic potential. 
The quadratic form defining the potential has two double eigenvalues. To see this we only have 
to calculate explicitly the potential (a, Ad 5 (r)) in the case when g G SU(2). The elements of 
SU(2) are matrices of the form 

9=( 91 + 2 92 93 + m ), det ( 5 ) = W = 1. 

1=1 

If we take 

( ia b + ic\ , ( i \ 

T= {-b + iC -ia) and (7= V -ij' 

we get 

2(a,Ad g (r)) = -Tr (a^" 1 ) 

= 2a(^ + ff2 - ff - ffl) + ^b(-gig 4 + 929z) + M9i93 + 929a)- 
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The matrix of the above quadratic form has indeed two double eigenvalues A = 2||r|| and p = — A. 
Thus the system (T*G,u c , H), where the Hamiltonian is given by (jSJ, can be considered as 
a generalized C. Neumann oscillator. Its configuration space is the Lie group G instead of the 
sphere. 

The fact that the quadratic form of the potential of (T*SU(2), u c , H) has two double eigen- 
values means that the system has two circular symmetries. Moreover, the system (T*G,lo c , H), 
where G is compact semi-simple, has generically two toroidal symmetries. Inspection of the 
Hamiltonian (jHJ) immediately shows that our system is preserved under the right action pl(g) = 
g ■ t of the maximal torus T T = expt,- and under the left action p\(g) = t ■ g of the torus 
T a = expt CT . Here t,- and t a denote the maximal toroidal subalgebras in g which contain the 
elements r and a, respectively. We have taken into the account the fact that the canonical 
symplectic form u c on T*G is left and right-invariant. The symmetries p r and p l will enable us 
to construct in Section |I] the symmetries of the generalized Maxwell-Bloch equation. 

Now we shall return to the partial differential equation (jHJ). We shall concentrate on the 
spatially periodic case, which means that we shall stipulate 

g(t, x + 2ir) = g(t, x), r(x + 2ir) = t(x). 

In a somewhat more graphic way, the equation © can be written as 

{gt.g^ 1 )^^) = --(gtg^i^x - e) - g t g~ l {t,x + e))\ e ^ Q + [a, Ad g(t>x) (r(x))]. 

Consider the function g(t,xo) at a fixed value xq of the spatial coordinate. We can think 
of g(t, xq) as of the position of the generalized C. Neumann oscillator moving in G. Its ac- 
celeration depends on the potential [a, Ad 9 (r)] and on the velocities (gtg" 1 )^, %o ± e) of the 
neighbouring oscillators. It is thus natural to try to understand the equation (jEJ) as the equation 
of motion for a continuous chain of the generalized C. Neumann oscillators. The oscillators 
interact among themselves in a magnetic way. By this we mean that the acceleration of every 
given oscillator depends on the velocity and not on the position of its neighbours. 

Let us now construct the Hamiltonian structure of the generalized Maxwell-Bloch equation 
which will correspond to the above interpretation. The obvious candidate for the configuration 
space of the continuous chain is the set of maps g(x) : S 1 — > G, or, in other words, the loop 
group LG. The phase space is therefore the cotangent bundle T*LG. The natural choice for the 
Hamiltonian is the total energy of all the oscillators: 

H mb (g(x),p g (x)) = J Q||p 3 (x)|| 2 + (a,Ad g{x) (T(x)))\ dx. (12) 

But the Hamiltonian system (T*LG, to c , H m b) does not correspond to the generalized Maxwell- 
Bloch equation ©. It is easily seen that the equation of motion of this system is simply the 
equation (g t g)t(t,x) = [a, Ad g ( t)X )(r(x))] which describes the system of decoupled generalized 
C. Neumann oscillators. Therefore we have to modify this system in such a way that the 
magnetic-type interactions will be taken into account. We shall achieve this by perturbing the 
canonical symplectic form tu c by an additional term. 

Let be the cocycle on the loop algebra Lg defining the central extension Lq = Lq x R. 
Recall that f2 is given by 

n(Z(x),T](x)) = - I (£'(x),r)(x)), £(x),r)(x) € Lg. 
Js 1 

Let us denote by ui m the right-invariant 2-form on the loop group LG whose value at the identity 
is equal to O, that is (cD m ) e = f2. Finally let ir: T*LG — > LG be the natural projection and 
let the 2-form to m on T*LG be given as the pull-back LJ m = TT*(u} m ). We have the following 
theorem. 
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Theorem 1. Let the Hamiltonian H mb in the Hamiltonian system (T*LG,u c + au m , H mb ) be 
given by and let uj m be of the form described above. Then the equation of motion of this 

Hamiltonian system is the generalized Maxwell-Bloch equation (JBJ). 

Proof. The proof will be a modification of the proof of Proposition ^ Note first that the 
Killing form on q induces an Ad-invariant inner product on Lq given by the formula 

«£(*), q(s)»= / 77(aO)ds. 
Js 1 

By the symbol ((— , — )) we shall also denote the pairing between the elements of Lq and those 
of Lq*, as well as the induced inner product on the dual Lq*. With this notation the directional 
derivative of the Hamiltonian H mb in the direction (5 b , 5 c t) is given by 

((&H mb ,(5 b ,5 ct ))) = -{{[a,Ad g (r)f,S b )) + «<Wj», 

where the maps ft and b are defined in the same way as before. 

The formula jUJ) from the Proposition ^ is valid for every Lie group, hence also for the loop 
group LG. Therefore we get the following expression for the form oj m + cu> m : 

(w c + cu> m ) {g ^ g) ((X b , X ct ),(Y b , Yet)) 

= -((X ct , Y b )) + {{Ya, X b )) + ((p g , [X b , Y b })) - c(((X b ) x ,Y b }). 

The above two formulae now give us the following relation for the Hamiltonian vector field 
x H mb = (X b ,X ct ): 

(u c + ao m )( gyPg )({X b , X ct ), (5 b ,8ct)) 

= -((Xct, Sb)) + ((Set, X b )) + ((p g , [X b , 5 b })) - c(((X b ) x , 5 b )) 
= ((- X ct - c(X b ) x - {X b , Pg }, 6 b )) + ((S ct , X b )). 

Since the vectors 6 b and 6 c t are linearly independent, we get two equations for the components 
of the Hamiltonian vector field: 

X b = p\, X ct + c(X b % = [a, Ad 9 (r)] b . (13) 

Let ~f(t,x) = (g(t,x),p g (t,x)): R — > T*LG be an integral curve of the field Xfj mh given in the 
right trivialization T*LG = LGx Lq. Then we have 7 = (gtg~ l , (Pg)t) an d thus from ()13|) finally 

{9t9~ l ) t + c(g t g- l ) x = [a,Ad g (r)]. ■ 

4 A family of conservation laws 

In this section we shall construct a family of conservation laws for the generalized Maxwell-Bloch 
system. This family will correspond to a certain group of symmetries of the Maxwell-Bloch 
system. The symmetries form a loop group. Hence the corresponding family of conservation 
laws will consist of infinitely many functionally independent elements. 
Let, for the sake of simplicity, the loop t(x) in the equation 

{gt9~ 1 ) t + c(g t g- 1 ) x = k,Ad 9 (r)] 

be constant, r(x) = r. Denote by tr the maximal toroidal subalgebra in q which contains the 
element r. Let T T = expt,- be the corresponding maximal torus in G. We notice immediately 
that our equation is invariant with respect to the right action 

Pr(x){g{t,x)) = g(t,x) ■ r(x), r(x): S 1 -> T T 
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of the loop group LT T C LG. Our conserved quantities will correspond to the LT T - symmetry in 
the sense of Noether's theorem. 

Equivalently, the generalized Maxwell-Bloch Hamiltonian system (T*LG,lo c + ao m , H m b) is 
invariant with respect to the canonical lifting g of the action p onto the cotangent bundle T*LG. 
It is easily checked that, in the trivialization of T*LG by the right translation, the action 
g: LT T x T*LG -» T*LG is given by 

0r(x){g{x),p g (x)) = (g{x) ■ r(x), p g (xj) . 

This is an action by symplectomorphisms. Indeed, the lifting to the cotangent bundle of an action 
on the base space is always symplectic. But this right action is symplectic also with respect 
to the form u m , due to the fact that u m is right-invariant. Thus the action g is Hamiltonian, 
since T*LG is a simply connected space. Therefore this action has the moment map. Our 
conservation laws will be the components of the moment map corresponding to the action g. 
We shall prove the following theorem. 

Theorem 2. Let (T*LG,lo c + cw m ,H m f,) be the generalized Maxwell-Bloch system. Let £ 
Ltr be an arbitrary loop inir. Then the function : T*LG — ► R, which in the right trivialization 
of T*LG is defined by the formula 

Ft(g(x), Pg (x)) = J i (p g (Ad g (0) - fag' 1 , Ad g (C))) Ax, 

is a first integral of this system. In other words, let 
g{t,x): RxS 1 ^G 

be a solution of the Maxwell-Bloch equation ©. Then for every loop £(x) : S 1 — > t T , the quantity 

G c (g(t,x))= [ (gtg' 1 -^g x g-\ Ad g (0)dx 
Js 1 1 

is constant with respect to time t. 

Proof. The second part of the theorem follows immediately from the first part. Let 7(t) = 
(g(t),p g (t)): R — ► T*LG be a solution of the Hamiltonian system (T*LG,u c + ao m , H m b), and 
let g(t,x): R x S 1 — > G be the map given by the curve g(t) = 7r(7(t)): R — ► LG, where 
7r: T*LG — > LG is the natural projection. Then 

g{t,x): RxS 1 ^G 

is a solution of the Maxwell-Bloch equation @. Again we shall use the trivialization T*LG = 
LG x g* by the right translations. Denote by 

pr 2 : T*LG = LG xg*^ Lg* 
the projection on the second component. Then the Killing form on G yields the identification 

W2{l(t)) =p g {t) = (g t (t,x)g~ 1 {t,x),-). 
This clearly shows that 

*e(7(*)) = / (gig" 1 - %9*9~\ Ad g (e))dx = G^{ g (t,x)) 

Js 1 1 

for every t. 
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We shall now show that are integrals of the Maxwell-Bloch Hamiltonian system. We have 
already mentioned that our integrals are the components of the moment map corresponding to 
the action g of TL T on the Maxwell-Bloch system. 

Let £(x) € Ltr be an arbitrary element in the Lie algebra Ltr- Denote by £ the infinitesimal 
action of £ on T*LG. By definition, the ^-component of the moment map fi: T*LG — > (Ltr)* 
is given by the formula 

Kd,Pg)(0 = H^(g,p g ), 

where : T*LG — > R is the function whose Hamiltonian vector field is the infinitesimal action £. 
This means that 

(&H$)(g,p g )(tib,8ct) = (Uc + CUJm)(£, (6b,5 ct )) = UJ C (£, (5 b ,5 ct )) + CU m (£, (6b, Set)) 

for every tangent vector (5b, 5a) £ T^ g ^ Pg ^(T* LG). Let the functions H^,H™: T*LG — ► R be 
given by the formulae 

(dfl"|)( s ,p g )(4,^ct) = (u c )(g, Pg )(£, (5 b ,5 ct )), (AH^) { g )Pg) (5 b ,5 c t) = (u m )(g,p g )(€, (5b,5 c t))- 

Then 



H ( = Hl + cH^. (14) 

Here we ignore the non-relevant indeterminate additive constant, and we shall continue to do 
so below. For the canonical form we have lv c = d9, where 6 is the tautological 1-form on T*LG. 
Thus we have 

( dH l)(g,p g )( S b,6 ct ) = (d(9) (9iP9 )(£, (5 b ,5 ct )), therefore H^(g,p g ) = 0( g , Pg )(€}- 
The expression of the infinitesimal action in the right trivialization is given by 

1(9, Pg) = ^ s _ o Qh(s)(g,Pg) = (Ad g (£),pg). 

Here h(s): (— e, e) — ► LT T is a path, such that h(0) = e, and ^| s=0 ^( s ) = C £ Lt-. From the 
definition of the tautological 1-form we now have 

Ht(g,Pg)=p g (l)= f Pg (Ad g (0)dx, (15) 
Js 1 

which gives the first summand in (|14j). 

We claim the the magnetic component H™ is given by the formula 

H?(g,Pg) = -hg x g-\ Ad g (0)) = - \ [ (g x g~\ Ad g (0)dx. (16) 
^ * Js 1 

We will calculate the derivative of H™ at (g,p g ) in the direction (5b, 5a) E T^ g . Pg ^(T*LG). Let 

(9(s),P g (s)) be a path, such that (c/(0),p g (0)) = (g,p g ) and aj| s=0 (ff(s),P s (s)) = (5b,5 c t) in 
the right trivialization. Derivation gives 

<WW)0M<*) = £ _„/ <0(*)*0OO~\Ad, (s) (O) 

US s— U J 51 



(((5^, Ad fl (0) + <fe<T\ [<5 6 , Ad fl (0]» dx 
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(S b , {Ad g (0)J ~ (fib, fag' 1 , Ad g (C)])) dx 
2 / (5 b ,[g x g-\Ad g (0})dx. 



On the other hand we have 



( w ™)(s,p,)(£> (h, Set)) = / ((Adg(£))x,5fe)dx = / (5 b ,[g x g 1 , Ad s (£)]) dx, 
that proves (|T5|) . 

If we now put together the formulae (|14|). (|15|) and l|16j) . we finally get 



Hi(g, Pg ) = j (p g (Ad g (0) - \(g x g' x , Ad 3 (£)>) dx 



that proves the theorem. ■ 

Our conservation laws can be expressed in a more compact and suggestive form. Let us 
choose an arbitrary point xo G S 1 and let ip n (x): S 1 — > R be a sequence of functions which 
approximates the Dirac delta function 5 Xo on S 1 . For every function / : S 1 — > R we then have 

lim / /(x) • <^ n (x)dx = / /(x) • o" Zo (x)dx = f(x ). 

n-*ooJ sl J sl 

Let £ G t T be an arbitrary element and let £ n (x) G Lt T be elements in the loop algebra defined 
by £ n (x) = ip n ( x ) ■ £■ Then it follows from the above theorem that for every n G N we have 

J^igtg- 1 - ^g x g~ l , Adg(£ n )) dx = J ([g t g~ l - ^g x g~ 1 ,Ad g (i))(x) ■ <p n (xf) dx = A n . 

The sequence {^4 n } ne N of constants is clearly convergent and its limit is 

(gtg^ 1 - ^g x g~ l , Ad fl (|))(x ) = y ^ {{gig" 1 - |ff a .ff" 1 ,Ad s (|))(x) • dx. 

This proves the following corollary. 

Corollary 1. For every element £ in i/ie maximal Abelian subalgebra ty- C 5 and for every point 
Xo € 5 , £/ie quantity 

Xx (g(t,x)) = {gig" 1 - ^ft^ -1 , Adg(£))(x ) 

is constant along every solution git, x) : R x 5 1 — ► G 0/ i/ie generalized Maxwell-Bloch equation. 
From the Ad-invariance of the Killing form (—, — ) we see that for every £ G £/ie quantity 

Xx = (g~ x gt - ^g~ 1 g x ,i)(xo) 

is constant. To put it more briefly, let p^ : g — > t r oe i/ie orthogonal projection. Then for every 
solution g(t, x) : R x S 11 — > G and /or every point xq G 5 1 i/ie projection 

Ptr ((f^ft - ^g^gx) (x )) 

is a constant element in tr. 
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Consider again the Hamiltonian function H m t,(g,p g ) of the Maxwell-Bloch dynamical system 
(T*LG,lu c + cw m , H m b). We see that the left action p l s (x)(g(x)) = s(x) • g(x) of the subgroup 
LT a C LG preserves the Hamiltonian H m b- Here T a is the maximal torus in G whose Lie 
algebra t a contains the element a. But the group LT a does not preserve the magnetic part uj m 
of the symplectic structure oj c + ao m . Recall that the magnetic part is right invariant, but it is 
not left invariant - not even with respect to the action of LT a . It is easily seen that only the 
torus T a C LT a (containing the constant loops) preserves the symplectic structure, and hence 
the whole Hamiltonian system. Calculations, similar to those above, only simpler, give us the 
proof of the following proposition. 

Proposition 2. Let r] € t a be an arbitrary element. Then for every solution g(t, x): K x 5 1 — > G 
the quantity 



of t a is constant along every solution. Here pt : Q — > t CT is again the orthogonal projection. 

It is well known that the Maxwell-Bloch equations are integrable. This is also true for the 
generalized Maxwell-Bloch equations. In particular, they satisfy the zero-curvature condition 



for the Lax pair 

U = -{- zo + g t g~ l ) and V = -za + gtg" 1 - -Ad s (r). 

The proof is a matter of trivial checking. We have to note that the conserved quantities con- 
structed above do not comprise a complete system of first integrals. But they can be used to 
reduce the system. An interested reader can find various results concerning the integrability of 
the Maxwell-Bloch equations in the references |121 1131 1141 115j etc. 

As we mentioned above, the conserved quantities associated with the action of LT T can be 
used to construct symplectic reductions of the system (T*LG, us c + cw m , H m f,). One can see that 
the nature of the symplectic reductions p~ l {a)/LT T depends quite heavily on the choice of the 
level a G Ltr of the moment map. We intend to study some cases in another paper. 

For the case of the original Maxwell-Bloch equation, these symplectic quotients can be 
thought of as infinite-dimensional analogues of the following situation. Let (T* SU(2),uj c , H), 



be the Hamiltonian system of the C. Neumann oscillator on S s with two circular symmetries, or 
equivalently, the Maxwell-Bloch system for the solutions which are constant with respect to the 
spatial variable. Then the right action of the maximal torus T T = exp t,- on our system is Hamil- 
tonian. Let fi: T*G — > t,- be the corresponding moment map. The quotient (/i _1 (a)/T T , u sq , H sq ) 
is the system (T*S 2 ,lo c + aw m , H sq ), where the induced Hamiltonian is given by 




is a constant. This means that the element 




V t -U x + [U,V} = 



H(g,p g ) = ^\\p g \\ 2 + (a,Ad 9 (r)) 



H sq (q,P q ) = ^\\p g \\ 2 + (<r,q) 
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and q = Ad g (r) S O t = S 2 C su(2). The magnetic perturbation io m of the symplectic form is 
given by 

(u m )( q>Pq )(X 1 ,X 2 ) =p q ([Kx,K 2 ]), Xi = [Ki,q] eT q O T = T q S 2 . 

The reduced system describes the charged spherical pendulum which moves in the magnetic 
field of a Dirac monopole situated in the centre of the sphere S 2 . More details and the proof 
can be found for example in |16U17j . The quotient construction is easily generalized to the cases 
when SU (2) is replaced by a compact semi-simple G. The sphere is then replaced by the general 
co- adjoint orbit O r of r and the magnetic form uj m is just the well-known Kostant-Kirillov 
symplectic form which is a part of the Kahler structure on O t . 

5 Concluding remarks 

We have described a new Hamiltonian structure of the Maxwell-Bloch equations. The key step 
in the construction is the observation that the Maxwell-Bloch equations can be considered as 
the equation of motion of a continuous chain of C. Neumann oscillators on the three-dimensional 
sphere. These oscillators interact through forces of the magnetic type. Our Hamiltonian struc- 
ture is derived from the well-known Hamiltonian structure of the C. Neumann oscillator. What 
has to be added in the case of the Maxwell-Bloch equations is the component which accounts 
for the magnetic-type interactions among the oscillators in the chain. This is achieved by means 
of the perturbation of the canonical symplectic structure by a topologically non-trivial magnetic 
term. This term is studied in more detail in [Z|. Among other things, we show in [2j that in the 
Lagrangian formulation this term gives rise to essentially the same topologically non-trivial term 
as the one appearing in the Wess-Zumino-Witten theory. In this paper we have concentrated on 
finding a family of conservation laws of the Maxwell-Bloch equations. The task was facilitated 
exactly by the use of our particular Hamiltonian structure. 

It is important to note that the Maxwell-Bloch equations might very well be endowed with 
additional Hamiltonian structures, different form the one described here. One such structure is 
given by Holm and Kovacic in |18| . It can be easily seen that the two structures are inequivalent. 
The symplectic structure presented in ^H] does not include derivatives with respect to the 
spatial variable, while ours does. The existence of these two Hamiltonian structures suggests 
the possibility of the existence of a bi-Hamiltonian structure on the Maxwell-Bloch system. 
We intend to address this issue in another paper. Another interesting avenue to the study of 
Hamiltonian structures was suggested to the author by one of the referees. He/she points out 
the important fact that the reduced Maxwell-Bloch equations can re rewritten in the form of 
a nonlinear von Neumann equation. This is shown in [T^]. A treatment of the nonlinear von 
Neumann-type equations and of their importance to quantum physics can be found in [20]. Some 
of the nonlinear von Neumann equations indeed have the Hamiltonian formulations of different 
sorts. An interesting example pointed out by the referee is the generalized Hamiltonian structure 
of Nambu, described in . 

Let us briefly outline a somewhat more traditional way of describing the Hamiltonian struc- 
ture of the reduced Maxwell-Bloch equations. The reduced system equations considered in |19j 
(and in other sources) is the system of three ordinary differential equations 

ii\ = — Au2, A = const, 

ii 2 = Atii + k£uz, k = const, 

u 3 = -k£u 2 , £ = (2/fcr)sech(Vr). (17) 

Since the variables u\, u 2 and u 3 are related by the constraint u\ + u 2 , + u 2 = 1, the phase 
space of this system is the two-dimensional sphere S 2 . We can view S 2 as a co-adjoint orbit of 
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a non-zero element in the Lie algebra su(2). If we introduce su(2)-valued maps g: S 2 — > su(2) 
and 0: S 2 x R -> su(2) 

IU3 ui+iu 2 \ O(t) — I ( ^ n£(t) 

-U1 + 1U2 —iuz J ' 2 \—K,£(t) —iA 

then the system (|17j) can be written in the form 

g = [n,g\. (18) 

Let us recall the Kostant-Kirillov symplectic form on the arbitrary co-adjoint orbit O C Q given 
by the formula 

(u kk ) m (X,Y) = m([a,b]), X = [m,a], Y = [m, b] E T m O. 

In the case of g = su(2), we can identify su(2) and the dual space su(2)* via the Killing form 
(a, (3) = —\ Tr(a • (3). Thus, on S 2 C su(2) we get the expression 



{u kk ) e {X,Y) = (g,[a,b]) = --Tr(g-[a,b]), X = [g,a], Y = [g, b] G T e S 2 . 



Let us now consider the function H : S 2 x R — > R given by the formula 

H(g,t) = (g,n(t)). (19) 

A short calculation shows that the reduced Maxwell-Bloch equations (|18|) are the equation of 
motion for the Hamiltonian system (T* S 2 ,oj kk , H) with the time-dependent Hamiltonian given 
by the formula (|19[) . 

We conclude this remark by observing that, in one way or another, the reduced Maxwell- 
Bloch equations lead to Hamiltonian structures which are different from the one considered in 
our paper. The main reason for this probably lies in the fact that in the case of the reduced 
equations, the electric field £ is not an unknown, but it is assumed to be a known function of 
time, for example £(t) = (2/kr) sech(t/r). 

In this paper, we studied the Maxwell-Bloch equations in their slowly varying envelopes 
approximation and in the sharp line limit, that is, without the inhomogeneous broadening. 
When the inhomogeneous broadening is taken into account, the resulting equations preserve the 
families of symmetries described in this paper. These symmetries can then be used for reduction 
and the resulting reduced equations would be analogous to those described above. However, the 
equations with broadening do not have a Hamiltonian structure directly analogous to the one 
described in this paper. For any Hamiltonian structure of the equations with broadening, the 
symplectic structure would have to be substantially different from ours. At this moment the 
author does not know, whether it would also carry over to the case without the assumption of 
the slowly varying envelopes. 
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